Conductivity Due to Classical Phase Fluctuations in a Model For High-Tc 

Superconductors 
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We consider the real part of the conductivity, cr\-aa{ijj), arising from classical phase fluctuations in 
a model for high-Tc superconductors. We show that ai-aaduj ^ below the superconducting 
transition temperature Tc, provided there is some quenched disorder in the system. Furthermore, 
for a fixed amount of quenched disorder, this integral at low temperatures is proportional to the 
zero-temperature superfluid density, in agreement with experiment. We calculate ai-aa{oj) explicitly 
for a model of overdamped phase fluctuations. 
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Because of their high transition temperatures, small 
coherence lengths, and low superfluid densities, the 
cuprate superconductors are strikingly influenced by 
phase fluctuations of the superconducting order param- 
eter. Such flufijtuations are largely responsible for flux 
lattice meltingtl and vortex glassci transitions in a finite 
magnetic field. Jn addition, they strongly affect the zero- 
field transitional, and possibly also the superconductinc 
transition temperature itself in underdoped materialsQ. 
Phase fluctuations also influence the transport proper- 
ties of the high-Tc materials. For example, the finite- 
frequeriGy conductivity shows a fluctuation-induced peak 
near tJB. 

Recent measurements in the most anisotropic high-Tc 
materials, such as Bi2Sr2Ca2Cu208+5, (BSCCO) have 
found that the real part of the conductivity, ai{uj), not 
only has a peak ijcar Tc but also remains quite large 
even far below TJb. Published measurements are avail- 
able typically at frequencies in the 30-200 GHz rangeSfl. 
In some cases, ai{uj,T) at such low temperatures ex- 
ceeds the peak values observed near T^. These large 
values occur over a broad concentration range varying 
from overdoped to optimally doped to underdoped. Fi- 
nally, large low-temperature values of ai{uj,T) are cor- 
related with a relatively large in-plane superfluid density 
ns{0) at T = 0. 

In this paper, we show that if the phases are assumed 
to fluctuate classically, then the presence of quenched 
disorder will produce a low-temperature background in 
CTi(w) which has many similarities to that observed in 
experimcnijj. Specifically, we argue that in a sample 
with quenched disorder there must inevitably be a low- 
temperature background ailujjT) arising from classical 
phase fluctuations, whose frequency integral is related 
to the zero-temperature superfluid density, ns{0). While 
the exact frequency dependence of this background de- 
pends on the particular dynamics obeyed by the phase 
fluctuations, the frequency integral is independent of the 
dynamics, but depends only on the assumption that the 
phases fluctuate classically. These fluctuations thus pro- 



vide an alternative possible source of conductivity back- 
ground, in addition to the gapless quasiparticles which 
should exist in a d-wave superconductorLl. 

In order to model the phase fluctuations, we adopt a 
classical XY model on a tetragonal lattice with different 
couplings in the ab and c directionaid. We consider a 
lattice model of a superconductor, such that each lattice 
point i is characterized by a phase 9i. In the presence of 
a vector potential, the interaction Hamiltonian for this 
system is given by 



H 



= - ^ Jij-W cos{e.i - Oj + Aij) 
^ Jkt,± cos(e'fe -9e + Am), 



(1) 



where the first sum runs over nearest neighbors in the 
ah plane and the second over bonds in the c direction; 
(T^j 1 1 and Jijjj_ are the couplings between lattice points 
in the ab and c directions. The gauge factor is Aij ~ 

(27r/$o) // -A- • d^, where A is the vector potential, $0 = 
hc/q is the flux quantum, and q is the magnitude of the 
charge of a Cooper pair. 

The diagonal components Us-aa of the superfluid den- 
sity tensor are given by ris-^aa — 
[{m*c'^/{Vq'^)]{d'^F/dA'^)A^=o- Here m* is the mass of a 
Cooper pair, F the Helmholtz free energy, V the sample 
volume, c the speed of light, and Aa is a flctitious uniform 
vector potential applied in the a direction in the presence 
of periodic boundary conditions in all three directions. 
Using F = -kBTlnUfL^J^"' dOiexpi-H/kBT), where 
N is the number of lattice points, we obtain 



-(7l;c 



■ 72;, 



(2) 



Here ji-a — [{2Tr/^o)aa]'^Ea, where aa is the lattice 
constant in direction a (ax — ay ^ a^ for the tetrag- 
onal lattice), <I>o = hc/q is the flux quantum, and 
Ea — (X](ij)||Q '^ij "^os %)j the sum running over bonds in 



1 



the direction a. Similarly, 72;a = [V^ / {c^kBT)]{Ja{tf) , 

where Ja{t) = {aa/V)J2{tj)\\ailJij/^)^^^^i3' ^ij ^^i- 

9j , the sums are carried out over distinct bonds in the a 
direction, and the triangular brackets denote a canonical 
thermal average. 

The key point is to note that J'a{t) is the volume- 
averaged supercurrent density in directioijua. From this 
connection, we can use the Kubo formulaHl, in the clas- 
sical limit, together with eq. (|^), to obtain an expression 
for (Ti;Q,a(w), the real part of the long-wavelength con- 
ductivity in the direction a: 



0'l;aQ(w) 



where L- 



V 

VkeT 
- qJij/h and 



dtCOs{ujt){Ja{t)JaiO)) 



dtcos(wt)(X„(0)Ia(t)), (3) 



ij)\\a 



Ic-ij sin 9ij{t). 



(4) 



Integrating eq. (||) over frequency leads to 
ai,aa{uj)dcj = [Vn/{2kBT)]{J^). Substituting this 
expression into eq. (||), and using CTctc((— w) = aaa{^), we 
finally obtain 



Sl-a = i 
JQ 



TTC 



/ 27rao 



Ea 



(5) 



We first show that, for an ordered system, S'5'.„(r) — 
at T = 0. Suppose that Jij = J a for all bonds in the 
a direction. Then the phases 9i are all equal and, with 
periodic boundary conditions, Ea — NJa- Similarly, in 
the limit T -> 0, id^F/dAl)A^=o = {d''H/dAl)A^=o = 
{2TTaa/^o)'^ N Ja, from direct evaluation of the second 
derivative, with periodic boundary conditions. Hence, 
both terms on the right-hand side of eq. (^) are equal, 
and S'?.„(0) = 0. 

We can use this result to obtain some analytical re- 
sults for S'i;q(T = 0) in the presence of quenched disor- 
der. Since the right-hand side of eq. (||) vanishes for the 
ordered case, <„„(0)gVw*c2 = {2naa/^o?E°{0)/V , 
where the superscript refers to the ordered system. Us- 
ing this equivalence, we may rewrite Si-a{T) as 



2m* 



(6) 



If the disordered system is unfrustrated (all > 0), all 
9i are still equal at T = 0, and Ea{T = 0) is controlled 
by the average bond strength X,: Ea{0) = NJa. In view 
of this fact, we choose reference system to have = Jq, 
so that Ea{0)/E°{0) = 1. 

To further analyze this regime, wCrHse a generalization 
of a theorem proved by KirkpatriclJl3, which maps the 



spin- wave stiffness constant of a random Heisenberg fer- 
romagnet onto the conductance of a disordered conduc- 
tance network, and generalized to the superfluid density 
of XY systen:is-|(the analog of the spin- wave stiffness con- 
stant) in Ref.t3. The generalized theorem, as applied to 
the present anisotropic geometry, states that 



jaoi. 



(7) 



Here gaa and gaa are components of the conductiv- 
ity tensor of fictitious random conductance networks in 
which the bond conductances are Jij and J°j, the bond 
strengths of the real and reference systems. 

Using this theorem, we can evaluate S'i;a(0) in an 
isotropic system for weak disorder. We retain our choice 
of the Ja for the reference system. Eor an isotropic 
conductance network in d dimensionaij, g/g — 1 — 
{{S9)'^)dts/[dg'^] + 0{{Sg)^), where 5g = g-g and {...)dis 
denotes an average over configurations of the quenched 
disorder. Hence, making use of eq. (j^), we obtain 

n.;aa(0) = <;„a(0) (l - ^'^'^ 'jfj , whcre d = 2 or 3 

is the dimensionality and J is the average bond strength 
of a bond in this network. 

Substituting this expression back into eq. (^) gives our 
final result for the integrated fluctuation conductivity in 
the case of weak disorder: 



fi;QC((a;, T — 0)duj 



7^g^7^s;Qa(0) ((J 



JY)dis 



2m* 



J^d 



(8) 



Here we have used the fact that, ns-aa{0) and n°.aa{Q) 
are equal to lowest order in ((J — J)'^)dis- The dimen- 
sions d = 2 or d = 3 would be appropriate for the lim- 
iting cases of no intcrplanar coupling (d = 2) or a fully 
isotropic system (d = 3); presumably a material such as 
Bi2Sr2Ca2Cu208+5W0uld behave similarly to the d ^ 2 
case. 

The resul|fc|p (|^) has several similarities 
to experimentsffl. Most strikingly, at fixed magnitude 
of the disorder parameter {{J — J)^)dis/J^, Si-a{T — 0) 
is predicted to scale with the low-temperature superfluid 
density, in agreement with experiment. Of course, the 
experiment is carried out at a specific frequency, while 
the relationship (||) applies to a frequency integral, but 
presumably CFi-aa{'-^) behave similarly for any given fre- 
quency. Note also that the result does not require the 
presenco-or absence of short range order in the bond 
strengthlla, but only that the bond distribution of bonds 
be macroscopically isotropic in d dimensions. 

It is of interest to make a numerical estimate 
of the integral (^) for parameters appropriate for 
Bi2Sr2Ca2Cu2 08+5. From the London equations, we 
may write irq^ris-aa/ {"^m*) = c^/(8A^). If we take the 
in-plane penetration depth as 2000A and we use a low- 
temperature conductivity ai{uj,0) of about 10^ fi^^m^^ 



2 



at u! 140 GHz, as suggested in Ref.El, and we further as- 
sume that this value remains roughly constant to low fre- 
quencies and cuts off at only slightly higher frequencies, 
then relation could be satisfied for (( J — J)^)/ ~ 
10" ^ This is at best an order-of-magnitude estimate, 
but it suggests at least that the relation (||) is not ruled 
out by experiment, since mean-square fluctuations in this 
range seem reasonable. 

The present results are entirely independent of the spe- 
cific dynamics: any classical dynamics will give the same 
value for Si-^a{T). Nevertheless, in order to illustrate 
a possible behavior for ai{io), we have approximated the 
dynamical response of the phase by the equations describ- 
ing a coupled array of mierdamped Josephson junctions 
[see, for example, Ref. (cj)]. In this picture, each lattice 
point is connected by an overdamped Josephson junction 
which carries three current contributions in parallel: a su- 
percurrent Ic-^ij sin^y, a normal current through a shunt 
resistance Rij, and a Langevin noise current lL;ij{t) rep- 
resenting the effects of thermal fluctuations. 

We will consider the anisotropic limit, in which the 
bond strengths vanish for ii\\c. In this limit, the 
supercurrents from bonds in different layers are uncorre- 
lated and eq. (^) reduces to (for conductivity lTl■_xx{'-^^T) 
parallel to the a axis) 



Sigma ^ 

300 



1 



dtCOs{ujt){Xa{Q)'Ia{t)), (9) 



where Ns is the number of lattice points in one layer, and 
the sum deflning 1^ [eq. (|^)] runs over lattice points in 
a single layer. We have evaluated this average by solv- 
ing the coupled pipsephson equations, using the method 
described in ref.Ej. 

Figs. 1 and 2 show ai-^xx{'^,T) for this model in two 
cases. In Fig. 1, all the critical currents and shunt re- 
sistance have the same values, Ic and R. In Fig. 2, the 
Ic;ij are uniformly and independently distributed on the 
interval (0, 2/^), but all the i?'s remain identical. In both 
cases, time is expressed in units of ?i/{qRIc), frequency 
in units of qRIc/h, current in units of Ic, ksT in units 
of hic/q, and therefore ai-xx{'-^,T) in units of l/{Raz). 

The results for the two cases are strikingly different. 
For the ordered lattice, there is a fluctuation peak in 
ai{uj,T), more prominent at smaller frequencies, cen- 
tered neat-the lattice phase-ordering temperature Tc ~ 
0.95hIc/<Jn. [The peak is probably shifted away from 
this value by flnite-size effects in our calculations.] For 
T < Tc, at all values of w, ai{(jj,T) falls off sharply to- 
wards a very small value at T = 0, consistent with the 
prediction that Si-aa{0) = 0. [The decreasing character 
of cri{uj, T) is most evident in the semilog plot.] 




FIG. 1. Fluctuation conductivity ai{ui,T), plotted as 
a function of temperature T for several values of the fre- 
quency u, calculated for an ordered 32 x 32 lattice with peri- 
odic boundary conditions. The inset shows same results in a 
semilog scale, demonstrating that 0-1(0;, T ^ 0) ^ 0. Units 
are as described in the text, and lines are guides to the eye. 

In the disordered lattice, ai{uj,T) still has a strong 
ffuctuation peak, but in addition has a broad background 
for all T < Tc, which rolls off at larger frequencies. The 
frequency-dependence is probably Lorentzian, as would 
be expected if the current-current correlation function in 
(^) decays exponentially in time. We expect that the 
relaxation time t(T) should be of order Ti/ {2eRIc)- The 
semilog plot shows clearly that ai{uj,T) remains flnite 
even as T ^ 0. 

The disorder-induced broad background in our calcu- 
lations may appear rather small (much weaker than the 
peak near Tc). But this background is calculated in units 
of the strongly temperature-dependent counling energy 
J. According to one model, for examplctJ, J oc A^, 
where A is the mean-fleld energy gap, a quantity which 
decreases to zero at the mean-field transition tempera- 
ture Tco . When the temperature-dependence of J is prop- 
erly included, the fiuctuation peak may well prove to be 
much smaller than the low-temperature background, es- 
pecially when the low-temperature is relatively large; 
this behavior would be in agreement with experiments. 
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cally and numerically, that in a superconductor with 
quenched disorder, the finite-frequency electrical conduc- 
tivity ai{uj,T) remains non-zero at low temperatures if 
the order parameter has phase fluctuations which can be 
treated classically. If this assumption is satisfied, this re- 
sult is quite general, and independent of the dynamical 
equations obeyed by these fluctuation. For weak disor- 
der, the frequency integral of this fluctuation conductiv- 
ity scales proportional to the low-temperature superfluid 
density, in agreement with recent microwave experiments 
in Bi2Sr2Ca2Cu208+5. 
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FIG. 2. (Ji((jj,r) plotted as a function of T for several 
frequencies in a disordered lattice (32 x 32 array), in which 
the critical current is uniformly and randomly distributed on 
the interval (0, 27c). The inset shows same results in a semilog 
scale. Lines are guides for eyes. 

We believe that the disorder-induced background con- 
ductivity originates as follows. In the absence of disor- 
der, phase fluctuations at low T will occur in the form of 
long-wavelength, low-frequency "phase phonons." These 
are loosely analogous to spin waves in ferromagnets, but 
may be underdamped or overdamped, depending on the 
dynamics. In the presence of quenched disorder, the den- 
sity of states of these phase phonons is affected by the 
scattering of these excitations off the disorder; in addi- 
tion, the modes of a given frequency, instead of having a 
unique wave vector, will form a wave packet with a spread 
in wave vector. Analogous behavior is found in the spin 
wave spectra of ferromagnets when there is quenched dis- 
order in, e. g., the exchange constants. We speculate that 
these damped modes produce the extra contribution to 
the ai{u;,T) at low temperatures. 

There are several issues which we have not consid- 
ered in the present work. For example, classical phase 
fluctuations might be expected to be frozen out by 
quantum effectS|— at low T, as may happen in low-Tc 
superconductorsim. But in some high-Tc materials, the 
existence of nodal quasiparticles may provide a normal 
background which would impede this quantum freeze- 
out. It remains an open question to what temperature 
classical phase fluctuations persist in such materials as 

Bi2Sr2Ca2Cu208+5- 

In conclusion, we have demonstrated, both analyti- 
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